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Abstract 

Some identities in law in terms of planar complex valued Ornstein- 
(^T) . Ulilenbeck processes (Zt = Xt +iYt,t > 0) including planar Brownian 

motion are established and shown to be equivalent to the well known 
Bougerol identity for linear Brownian motion (/?t,t > 0): for any fixed 



OO 



\Q . u > 0: 



■ \.(a \ ( law ') a 
Smh(P u ) = P(J o " ds exp(2/3 s ))) 



o , 

f^ | with (/3j,t > 0) a Brownian motion, independent of /3. 

These identities in law for 2-dimensional processes allow to study the 
distributions of hitting times T® = inf{£ : 6t = c}, (c > 0), T_ dc = 

K> , inf{£ : 9t ^ (— d, c)}, (c,d > 0) and more specifically of T_ cc = inf{£ : 

0t $ (— c,c)}, (c > 0) of the continuous winding processes 9t = 
Im(L -%*■), t > of complex valued Ornstein-Uhlenbeck processes. 
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1 Introduction 

The conformal invariance of planar Brownian motion has deep consequences 
as to the structure of its trajectories (see, e.g., Le Gall [21]). In particular, a 
number of articles have been devoted to the study of its continuous winding 
process (9t,t > 0): Spitzer [33], Williams [35], Durrett [16| . Messulam-Yor 
[2H], Pitman- Yor [SDJ, Le Gall-Yor [22J, Bertoin- Werner [5J, Yor [38], Pap-Yor 
[29J, Bentkus- Pap-Yor [4j. In this paper, we take up again the study of the 
first hitting times: 

T e _ djC = ini{t:9 t ^(-d,c)}, (c,d>0), 

this time in relation with Bougerol's well-known identity (see Bougerol [8], 
Alili-Dufresne-Yor [2] and Yor [39J): for fixed u > 0: 

sinh(/3 u ) = /% o u d scxp (2/3 s )) ) 

where ((3 t ,t > 0) is a Brownian motionzl, independent of (3. 
In particular, it turns out that: for fixed c > 0: 

# r< /3 = C a (c), (*) 

where J3 is a BM* independent of (6 U , u > 0), T$ = inf {t : (5 t = c } (C u t > 0) 

is a standard Cauchy process and a(c) = arg sinh(c) = log (c + y/l + c 2 ) , c G 

R. 

The identity (•) yields yet another proof of the celebrated Spitzer theorem: 

2 Q (law) n 
Vt > W , 



log t i->°o 

with the help of Williams' "pinching method" (see Williams [35J and Messulam- 
Yor [28]). 

Moreover, we study the distributions of T_ ^ c and T^ c c . In particular, 
we give explicit formulae for the density function of T_ and for the first 
moment of In {T < L C c ) . 



*When we simply write: Brownian motion, we always mean real-valued Brownian mo- 
tion, starting from 0. For 2-dimensional Brownian motion, we indicate planar or complex 
BM. 



The last section of the paper is devoted to developing similar results 
when planar Brownian motion is replaced by a complex valued Ornstein- 
Uhlenbeck process. We note that Bertoin- Werner [5] already made discus- 
sions of windings for planar Brownian motion using arguments related to 
Ornstein-Uhlenbeck processes. 

Firstly, we obtain some analogue of (*) when T@ is replaced by T c = 



nZ 



T_ cc = inf{£ : \6f\ = c}, the corresponding time for an Ornstein-Uhlenbeck 
process with parameter A. Secondly, we identify the distribution of T c . 
More specifically, we derive the asymptotics of E Tc for A large and for 
A small. 



2 The Brownian motion case 

2.1 A reminder on planar Brownian motion 

Let (Z t = Xt+iYt, t > 0) denote a standard planar Brownian motion, starting 
from xq + iO,XQ > 0, where (X t ,t > 0) and (Y t ,t > 0) are two independent 
linear Brownian motions, starting respectively from Xq and 0. 

As is well known (see e.g. Ito-McKean [20]), since Xq ^ 0, (Z t ,t > 0) 
does not visit a.s. the point but keeps winding around infinitely often. 
In particular, the continuous winding process 9 t = Im(J 4p-),t > is well 
defined. 

Furthermore, there is the skew product representation: 

ft A<7 

log \z t \ + ie t = -^ = (Pu + iiu) „ rf ,. , (i) 



J 8 



-Ht=R 



where (j3 u + «7 n , u > 0) is another planar Brownian motion starting from 
logxo + i0. For a study of the Bessel clock H, see Yor [36J. 
Rewriting ([1]) as: 

\og\Z t \ = p Ht ; 6 t = lHtJ (2) 

we easily obtain that the total a-fields a\\Z t \ , t > 0} and a{(3 u ,u > 0} are 
identical, whereas (j u ,u > 0) is independent from (\Z t \ ,t > 0). 



Figure 1: Exit times for a planar BM. This figure presents the exit times 
(a) T c e (ti doesn't matter because the angle is negative) and (b) T e _ dc for a 
planar BM starting from x + iO. 

A number of studies of the properties of the first hitting time (see Figure 

T e _ dfi = inf{£ : 9 t £ (-d, c)}, (c, d > 0), 

have been developed, going back to Spitzer |33| . 

In particular, it is well known (Spitzer |33| . Burkholder [9], Revuz-Yor |32] 

Ex. 2.21/page 196) that: 

E [(T e _ d j] < oo if and only if p < -^r- (3) 

Moreover, Spitzer's asymptotic theorem (see e.g. Spitzer [33j) states that: 

29 1 {law) (law) , . N 

' C x = 7 T o, (4) 



iOgt t-+oo 1 i 

where C\ is a standard Cauchy variable. 



2.2 On the Laplace transform of the distribution of the 
hitting time T e c = T^ c 

Now, we use the representation (J2J) to access the distribution of T c 9 (see Figure 
Ufa)). We define T c 7 = inf{t : j t ^ (— oo,c)} the hitting time associated to 



the Brownian motion (7t,£ > 0). Note that, from (J2J): 
H T e = T c 7 , hence: T d c = H~} , where 

H- 1 = inf {t :H t >u}= ds exp(2&) := A u . (5) 

Thus, we have obtained: 

Tt = A T2 , (6) 

where (A u , u > 0) and T c 7 are independent, since /3 and 7 are independent. 
We can write: f3 s = (logXo) + fis , with (/3s ,s > 0) a standard one- 
dimensional Brownian motion starting from 0. Then, we deduce from (jHj) 
that: 

T e c =xl(j\seMWi 0) )\ (7) 

From now on, for simplicity, we shall take Xq = 1, but this is really no 
restriction, as the dependency in Xq, which is exhibited in (171), is very simple. 
We shall also make use of Bougerol's identity [HI |2] and |3H] (p. 200), which 
is very useful to study the distribution of A u (e.g. |26l 127]). For any fixed 
u > 0: 

sinh(/3 n ) = /3a u = % "d S ex P (2/3 s )) ) ( 8 ) 

where on the right hand side, ((3 t ,t > 0) is a Brownian motion, independent 
of A u = J^dsexp(2/3 S ). 

Thus, from (jSJ) and (jSJ), and as is well known |32| . the law of (3 T i is the 
Cauchy law with parameter c, i.e., with density: 

^(2/) = TT] — 2T' 

7T(C Z + |/ Z J 

we deduce that: 

Proposition 2.1 For fixed c > 0, there is the following identity in law: 

smh(C c ) {l = ] P (T!) , (9) 

where, on the left hand side, (C c , c > 0) denotes a standard Cauchy process 
and on the right hand side, (f3 u , u > 0) is a one-dimensional BM, independent 
from T e c . 



We may now identify the densities of the variables found on both sides of 

flSD,i-e.: 

on the left hand side: 



on the right hand side: E 



/i c (argsinha;) = 

i ( x 2 

exp l 2W 



Vl+x- 



■h c (a(x)) ; 



where a(x) = argsinh(x). 
Thus, we have obtained the following: 



Proposition 2.2 The distribution ofTf? may be characterized by: 



E 



1 



V^f! 



exp 



x 



2T e 



1 



Vl + x tt(c 2 + log^^i + y/T+x)) 



x>0. 



(10) 



The proof of Proposition 12.21 follows from: a(y) = argsinh(y) = log(y + 
a/1 + y 2 ) and by making the change of variable y 2 = x. Let us now define 
the probability: 



Q c 



7r<r 
2Tf 



P 



The fact that Q c is a probability follows from ( liop by taking x = 0. Thus we 
obtain that c E[^ir /2Tf!] = 1, and we may write: 



Er 



exp 



X 



2T! 



y/TTx 1 + 4 log 2 (v/£ + y/l + x) 



Wx > 0, (11) 



which yields the Laplace transform of 1/Tj under Q c . 

Let us now take a look at what happens if we make c — > oo. If we denote by 
Tf = inf {t : (3 t = 1} the first hitting time of level 1 for a standard BM j3 and 
by A^ a standard Gaussian variable Af(0, 1), from equation (ITTj) . we obtain: 



lim E Qc 



-x/2T! 



E [e- xN2 ' 2 ) = E (e^/< 



(12) 



which means that : T 



e ( laW ) rp/3 



T[ . (At this point, one may wonder whether 

c— ^oo 

there is some kind of convergence in law involving (8 u ,u > 0), under Q c , as 
c — > oo, but, we shall not touch this point). 



From Proposition 12.21 we deduce the following: 



Corollary 2.3 Let <p(x) denote the Laplace transform < f7I)) . that is the Laplace 
transform ofl/2T® under Q c . Then, the Laplace transform o/l/2T c under 
P is: 



E 



cxp 



X 



2T! 



dw 



<p(w) 



(13) 



\/w — X 

Proof of Corollary 12.31 From Fubini's theorem, we deduce from (II ip that: 



E 



exp 



.r 



2T e 



which is formula (fl3"|) . 



y=xt 



v=l+t 



CO 



dy 



v 7 ^ 
°° dy 

00 dt 

" Vt 



I f x + y 

eXp t ~^" 



27rT„ e 

V9(x(l+t)) 







<il> 



i v«- 1 

tp{w 



ip(xv) 



V 



w — X 



a 



2.3 Some related identities in law 

This subsection is strongly related to [T5] . 

A slightly different look at the combination of Bougerol's identity (jHJ) and 
the skew-product representation ([1]) lead to the following striking identities 
in law: 

Proposition 2.4 Let (5 u ,u > 0) be a 1-dimensional Brownian motion in- 
dependent of the planar Brownian motion (Z u ,u > 0), starting from 1 + iO. 
Then, for any b > 0, the following identities in law hold: 



(i)H 2 



( la ^) rpP 



T 



a(b) 



/..\ a (law) n 
\\\) V T 5 = G a (6) 



,...s 5 (law) 

(111) V T s = |O a (b)|, 



where Ca is a Cauchy variable with parameter A and 8 U = sup s<u 8 s . 



Proof of Proposition 12.41 From the symmetry principle (see [3] for the 
original Note and [17] for a detailed discussion), Bougerol's identity may be 
equivalently stated as: 



sinh(/3 u ) = 5 AuW . 



(14) 



Consequently, the laws of the first hitting times of a fixed level b by the 
processes on each side of (|T4"|) are identical, that is: 



which is (i). 

(ii) follows from (i) since: 



P (law) 
1 a(b) ~ n T^ 



(law) 
Vu = 1H U , 



with (7 S , s > 0) a Brownian motion independent of (H u , u > 0) and (C u , u > 

0) may be represented as (j T 0,u > 0). 

(iii) follows from (ii), again with the help of the symmetry principle. 



□ 



Remark 2.5 Proposition ^. !& may also be derived from (iii) in Proposition 
\2.4\ Indeed, fore > 0, starting from the LHS of (iii), and letting N ~ A/"(0, 1) 
independent from Tf?: 



p [ e T 5 < c 



P{T & b <T°) = P{b<5 T e) 
P(b< ^T°\N\ 
b 



P 



/j^e 



<\N\ 



E 



dy e 



-y 2 /2 



b/y/Ti 



while, on the RHS of (iii): 



P(\C a{b) \<c)=2 



a(b) dy y =a(b)h 2 f c/a{b) dh 

71 



1 + h 2 ' 



(15) 



(16) 



/„ 7r(a 2 (6) + y 2 ) n Jo 

Taking derivatives in fTS\) and [ToV with respect to b and changing the vari 
ables b = \fx, we obtain Proposition ^. H 
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2.4 Recovering Spitzer's theorem 

The identity (ii) in Proposition 12.41 is reminiscent of Williams' remark (see 
1281). that: 



H TR = T logr , (17) 

where here R starts from 1 and S starts from (in fact, this is a consequence 
of (J2J) ). For a number of variants of ( IT7|) . see [37, 25J. This was D. Williams' 
starting point for a non-computational proof of Spitzer's result (j3j). We note 
that in (ii), T)f is independent of the process (9 u ,u > 0) while in (ITT)) T^ 
depends on (8 U , u > 0). Actually, we can mimic Williams' "pinching method" 
to derive Spitzer's theorem (J3J) from (ii) in Proposition 12.41 

Proposition 2.6 (A new proof of Spitzer's theorem) 

As t — > oo ; 6 T s — 9 t converges in law, which implies that: 

e T s - e t ) A o, (is) 



logt \ V* / t^oo 

and, in turn, implies Spitzer's theorem (see formula (Q) ): 

2 (law) 

t>t — > c-i. 



i'X- 



log t t- 
Proof of Proposition 12.61 From equation (ii) of Proposition 12.41 we note: 

1 (law) C a (b) (law) 

UrpS = — y Lj\. 



oo 



log b b log b b^ 

So, for b — yt we have: 

2 (law) „ 
o T s — y Ol- 

logt Vi b->oo 

On the other hand, following Williams' "pinching method", we note that: 

1 (q t& _ 9t ) <i=$ o, 



logt V v^ / t-*tx> 

since Z u = x$ + Z u and also, as we change variables u — tv and we use the 
scaling property, we obtain: 

•" \Jt Zu /'-><» Ui z°> 



Here, the limit variable is -in our opinion- of no other interest than its exis- 
tence which implies (1181) . hence (J3J). 

□ 

2.5 On the distributions of T e c = T^ c and T 6 _ cc 
Proposition 2.7 The asymptotic equivalence: 

(logt) P(T! > t) ^ (4c)/tt , (19) 

t4s a consequence, for r\ > 0, ^[(logXf)^] < oo if and only if r/ < 1 (where 
(■)_!_ denotes the positive part). 



Proof of Proposition 12.71 a) We rely upon the asymptotic distribution of 

ds 

Jza 2 



H t = f Q -^y which is given by [32]: 



-i% m Tf = M{t : ft = 1}, (20) 

(logrj z i->oo 



or equivalently: 

log £ (Joto) 
2v / Hy" t^oo 



N\, (21) 



where iV is a standard Gaussian variable A/"(0, 1). 

We note that, from the representation (J2J) of 6* i5 the result (1201) is equivalent 

to Spitzer's theorem [33] : 

20J (iau) (tow) , , 

1 — I — ► w = 1 T ^ \ 22 ) 

logt t-+oo J i 

where C\ is a standard Cauchy variable. 

b) We shall now use this, in order to deduce Proposition 12.71 We denote 

Sf = swp s<t 6 s = Sjj and we note that (from scaling): 

P{T e c >t) = P(Sl t < c) = P(0^ < c), (23) 
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since 7 and H are independent. Thus, we have (since Sj 
making the change of variable x = -%=)'■ 



(law) 



P{T» > t) 



2 -E 



cE 



\N\ and by 



c/vm y 

dx e~~ 




I V^ 6XP| 


f cWX\ 



(24) 



(25) 



(26) 



Thus, we now deduce from (|2"T|) that: 

*l±P{T* e >t)^J*cE[\N\] = Zc 

I \ 71 TV 

which is precisely (fT9|) . 

It is now elementary to deduce from (125]) that: for rj > 0: 

E[(\ogT 9 c )l} < 00 & 0<7?<1, 

since (I2"5j) is equivalent to: 

m P(logT c y > u) — > ( — ) • 

Consequently, Fubini's theorem yields: 

E l(\ogT e c )l] = dur] u^ 1 P(\ogT e c > u), 
Jo 

and from ( 126]) this is finite if and only if: 

du vP~ 2 < 00 ^77 < 1. 

So, E[(\og T e c )\] < 00 & < 77 < 1. 

D 

Now we give several examples of random times T : C(R + ,M) — > R + which 

may be studied quite similarly to T 6 e . 

For such times T, it will always be true that: Hj<(e) — T(j) is equivalent 



11 



to T(6) = At^), defined with respect to Z, issued from Xq 7^ 0. Using 
Bougerol's identity, we obtain: 



(law) 



sinh(/3 T(7) ) = (3 A = (5 {T{0)) . 



(27) 



where (J3 u ,u > 0) is a 1-dimensional Brownian motion independent of (/?, 7) 
(or equivalently, of Z). Consequently, denoting by h? the density of /3tM> 
we deduce from (127)) that: 



E 



x 



exp 



2nT(6) " 2T(6y 



VT+ 



xs 



h T {\og{x + Vl + x 2 )), 



f28) 



or equivalently, changing x in -y/x, we obtain: 



E 



exp( 



V^T^y 2T(0) 



vi + x 



fcr(log(>/x + VT+s)). (29) 



In a number of cases, Ht is known explicitly, for example: 

(i) 



T( 7 ) = T\ c & T(0) 



T , 



Gfs exp (2&) = 7* 



Thus: 



£ 



2^* 



exp(- 



.r 



2T 6 



d.c 



Vl + x 



h- d , c (log(Vx + y/T+x)), (30) 



where h-d c is the density of the variable j3 T ~t . The law of f3 T -r may be 

— d , c — d , c 

obtained from its characteristic function which is given by [32], page 73: 

A 2 



E 



exp(iA/3 T - 



d,c 



E 



exp{-—TZ djC ) 



cosh(|(c - d)) 
cosh(f(c + d))' 

In particular, for c = d, we recover the very classical formula: 

1 



£ 



exp(i\(3j 



cosh (Ac) 



12 



It is well known that [2"H |B] : 



E 



exp(i\(3 T i )) 



cosh (Ac) cosh(7rA £ ) 



-oo 
oo 



27r cosh( 



3 « A J/_ 



-oo 

OO I 



2vrcosh(ff 
1 



-dy 



oo 



2c cosh (- 



FT^- 



2c - 



Hence, the density of f3 T -t is: 



h. 



IX 



1 



2cy cosh(|p 



X7T X7T J 



and 



h_ CiC (log(v^+ V / T+i)) = (-) 



(31) 



c/ (v^ + v / l + s) c + (v 7 ^ + v / TTx) 



where C = f-- However using: 



{Vx + VT+x) c = {VTTx - Vx) c , 



(32) 



we obtain: 



fr_ c , c (log(v/i+ Vl + x)) = (-) 



So we deduce that (for c = d): 



(y/x + ^Jl + xY + {y/l + x - y/x)S ' 



(33) 



E 



1 x 

= eXp( -2T^ 



2ttT 6 



Cj \^l + xj {^/X + VI + XY + {y/l + X- y/x)< 



(34) 



□ 
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(ii) As a second example of a random time T, let us consider the time intro- 
duced in [34| . |1U| . exercise 6.2, p. 178 (we use a slightly different notation). 
Let (Pt,t > 0) be a real valued Brownian motion and define, for c > 0: 



T(0) = T° = inf it : sup 6 S - inf 6 S = c 

{ s<t s <* 

T{i) = T c 7 = inf -U : sup 7 S - inf 7 S = c 



s<* 



s<* 



Thus, from the skew-product representation ([TJ, U = 7^, by replacing 



u = T c , we obtain: 



Hp, = 77 =► 7? 



2;? 



<is exp (2/3 s ) = A^. 



Thus: 



£ 



1 , x 

exp 



2nT? 



2T e- 



v/TTx 



^(log^+VTT^)), (35) 



where h c is the density of the variable (3 T j. The law of /3 T7 may be obtained 
from its characteristic function which is given by P, ITU] : 



£ 



exp(zA/J, 



£ 



A 2 
exp(-yTJ) 



g V 2tt / 



(cosh(Af)) 2 (cosh(7rA^)) : 



.r 



-00 
00 



27rsinh(| 



-dx 



Ay A £H d 

27rsinh(^) c y 



c 2 sinh(- 



(36) 



So, the density of /5 T 7 is: 



M?/) 



2vry 



sinh( 



7TJ^ 
C ' 



47r y 

c e c — e c 
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and 



h c (log(y/x + VT+ 



where ( = -. Thus: 

3 r 



E 



x , 



4tt 



log(A/i+ y/1 +x) 



(y/x + VI + x) C ~ {y/x + VI + x)- 



exp( 



x 



2*7* 



2rr 



4tt 1 



log( V^ + VI + x) 



c 2 VTT^(V^+ VT+x)< - (VT+x - V^) c ' 



(37) 



We note that this study may be related to [31]; and more precisely P t j and 

T2 correspond to the variables C^ and C<i respectively (see e.g. Table 6 in p. 
312). 



□ 



Let us now return to the case of T e _ cc (example (i)). More specifically, we 
shall obtain its density function /(£). 



Proposition 2.8 The density function f ofT e _ cc is given by: 



f(t) 



V2< 



1 oo 



■* r( " t) l ,-iM, ., (I) 



fc=0 



r(2z/ fc ) Vt 



,"k i 



if 



(38) 



where M a ^(-) is the Whittaker function with parameters a,b. Equivalently: 



««) = ^E(-i)*V*^ 



fc=0 



i>* 



E 

n=0 



I> fc + n) 1/lV 



r(2z/ fc + n + 1) n! \2t 



(39) 



where v\. = j-(2k + 1 



Proof of Proposition 12.81 The following calculation relies upon a private 

note by A. Comtet [11] .We denote: 

(p C (x) = (y/x + V1 + X) C + (Vl + X - y/xf. 
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Noting: 



V 



we get: 



\/l + x = cosh - -<=>- y = 2 arg cosh(vT+"x) 



(f(;(x) = (sinh — + cosh 77) + (cosh — — sinh —p 

= 2 cosh — . 
2 



Thus, from (1541. we have: 



(40) 



II :=E 



x 



2vrT e 



: exp( 



2T 6 



■0 cosh I cosh H ' 



where ^ = 2c. However, expanding cosh 2jjL, we get: 



1 



_ Tvy_ 

e 2 * 



coshg 



1 wy 

1 + e * 



2E 

fc=0 



_ TV\j_ \ •• 7TJ/ 

-e * I e 2 * , 



(41) 



and from ( 14TJ) . we deduce that: 
77 



■0 cosh I 

4(-l) fc 



sinh I 



2-v/2 sinh | cosh | y cosh | 



2_ e -VkV 



jr 2 ^~!^. c~ T " i2 '' +j " / 

fc=0 
00 

= £ 

fc=0 

00 
f-f ipV2J2 sinh | cosh | 

fe=0 r V 2 2 

where u k = -^{2k+l). 

From (j40|) . we have 1 + x = cosh 2 | -<=>- x = sinh 2 |, thus: 



4(-l) fc 



tanh — e 



" --^3/ 



(tanh 



h 1 / 2 



' sinh I 



cosh I V a/1 + x 



.r 



1/2 



1 + x 



i/i 
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Moreover, we know that (see [1], equation 8.6.10, or 



7T 



1/2 



H 77— T-e-^ = g^_ 1/2 (coshy) 



2 sinh y 



Vk-l/V 



where {Ql{-)} is the family of Legendre functions and coshy = 2x + 1. So, 
we deduce: 



^-E^^^j^e^H-D- (-) 



4( ~°f i^ r x v ^ 

By using formula 7.621.9, page 864 in [18 



*w 2r(l + 2i/ fe ) e~ i7r ' 7 ' N ' 2 

u; ~ r(| + i/ fc + o Vs + 



^m^ T,?,V !t!) «W*). («) 



with: Z = |, z/fc = ^(2^ + 1), s = x + ~ and M v (-) denoting the Whittaker 
function, which is defined as: 

„„ , , b ,i i,„ r(26+l) ^r(i + 6-a + n) w ™ 
we have: 

(44) 
From ( 1421) and by changing the variable w = =i, we deduce: 

II = r ^(-1)' ^ + H r - exp f- W (x + -)1 M 1/2 Vh (w) 

f^VP ; r(2z/ fc + i)7 w P V 2 J J 1,2 ' VkK } 



k=0 



v^ / dt 2 , 1N , r(z/ fc + l) / l x\ % , ,1, 

(45) 
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By using the equations (I4ip and (I45p . we conclude: 



E 



1 



exp 



E 

k=0 

oo 

E 

fc=0 



00 dt 2 
dt 2 



.r 



2T e 



£ ^a/tt 



:-i)' 



-i) ; 



exp 



i> fc + 1) 
r(2z/ fc + i; 



T(^(2fc + !) + !) 

r(2£(2* + i) + i) 



■^--W(±' 

4t 2£/ *'" hV 2i' 



exp 



1 x \ , , ,1, 

4t 2£/ 2>4H( 2fe + 1 )^2t ; 



(46) 



Thus, the density function / of T e _ c c is given by: 



^_^ v ' tVo,,. li"i /I 5« l ^ ; 



72 
c 

v/2 



fc=0 
oo 



T(2u k + l)y/t 



2t 



D-if^g 11 ^"!,-^. 



1 



fc=0 
oo 

E< 

fc=0 



r(£(2fc + i) + iM 
. i/ fc r(i/ fc ) i _j, i, 

' ; — -r—pe 4* Ml „ ( — 

2u k T(2u k ) y/i ^ VkK 2V 



(47) 

t a i^ + i)\j t ) (48) 
(49) 



where the Whittaker function Mi .., (=r) is: 

2 )^fc ^ Zt ' 



Mi tf _ {2k+1) {-: 



1 ^ £(»+!)+£ __^ r ( * ( 2 fc + i) + i)- r(£(2fc + I ) + n) I f 1 

2-^ TYiL 



2t, 

1 \ v> " 

2f 



e 4 * ■ 



r(£(2* + i)) 



r(^(2A; + l) + l + n)n! \^2t 



=0 v 2c 



e 4t 



r(2z/ fc + 1) 
r(i/ fc ) 



n=0 

e-ir(2z/ fc )-^-^^ 



I> fc + n) 1 /l 



2t 



r(2z/ fc + 1 + n) n\ \2t 
T(v k + n) 



*>*) ^ (2^ + n)r(2z/ fc + n) n! \2t 



Thus, from (091) and (1501). we deduce (139 



(50) 
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1.5 2 2.5 3 3.5 4 4.5 5 

t 



Figure 2: Graph of fgg(t), with c = 2ir. 



a 

Next, we present the graphs of different approximations Jk,n{^) of /(£), 
in (139|) . where /k,jv denotes the sum in the series in (139|) of the terms for 
fc < A", and n < N. 

Remark 2.9 • Figure\^represents the approximation of the density func- 
tion f with respect to the time t (for K and N < 9), with c = 2n, 
whereas Figure represents the approximation of f with respect to the 
time t for several values of k and n, with c = 2tt. 

• From Figure [3J we may remark that the approximation K and N < 9 
is sufficiently good (comparing to the one for K and N < 100,). 

• For the case K and N < 9 it seems that locally, in a small area around 
0, f(t) < which is not right. This is due to the first negative (k = 1) 
term of the sum and due to the fact that we have omitted many terms. 
However, this is not a problem because it appears only locally. Similar 
irregularities have already been observed in previous articles JW$ p.275. 
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K = N = 9 

K = 0, N = 9 

K = 9, N = 

K=N=0 

K = N = 100 




Figure 3: Graph of fK,N{t) for several values of K and N, with c = 2ir. 

2.6 On the first moment of In (T_ cc ) 

This subsection is related to a result in |12| . 

Proposition 2.10 The first moment of ha. (T e cc ) has the following integral 
representation: 



E [In (T\ c )] = 2 ^ cogh In (sinh (cz)) + In (2) + c E , 



(51) 



where Ce = — r'(l) is the Euler-Mascheroni constant (also called Euler's 
constant). 



Proof of Proposition 12.101 Let us return to equations (J2]) and ([6]). So, for 
t = T® „, we have: 



— c,c 7 

Thus, for e > 0: 



77, 



T 9 ~~ ^""-c c 



T\.= An 



(52) 



^[(^ C ,c)1=^[(^ 
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Consider (S u ,u > 0) a Brownian motion, independent of A t . Then, Bougerol's 
identity and the scaling property yield (Q a denotes a gamma variable with 



. A7 .r> (law) _ < 

parameter a, and TV = 2yi/2) 



E[{sinh{B t )r] = E[(5 At ) 2£ ]=E[Al(5 1 f] 
= E[AI] E[(2g 1/2 ) £ ] 



because 



M(M1 



£ [At] (T 



e+i-l e 



r(i) 



r(l) 



rf.r 



Thus, for t = T_L „, we have: 



E 



sinh I By 



2c 



E 



A £ 



(2 e 



r(j) 



r(i) 



Recall that [SHE]: 



£ 



exp(zAi? T 7 )) 



£ 



A 2 
exp(-yT2 C)C )) 



cosh (Ac) 



and the density of (3 T t is: 



(53) 



h-c, e (y) = — 



2c J cosh(||) 
Thus, on the left hand side of ( | 53 |) . we have: 



y~K yir • 



E 



sinh I B* 



> T i 



2s 



2c cosh( 



— (sinh(?/)) 



■2 s 



2c- 



I f^km^^ 



OO -1 

rfz — T^ 

(l cosh( 



7T£ 

2 - 



(sinh(cz)) 



2s 
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where we have made the change of variable z = -. Hence, from f l53l) . by 
writing: 



E 
we deduce: 

r(|+f) 
r(i) 



^ 



E[(TtJ*\=E 



e £ ln ( T - c , c ) 



£ 



e e H^ce) 



d<3 



2 £ j cosh( 



7T£ 
2 - 



(sinh(cz)) 



2c 



and by removing 1 from both sides, we obtain: 



rQ + E) 
r(|) 



^ 



e £ H T -ccJ 



cte / (sinh(cz)) 



2c 



cosh( 



7TZ ' 

2 / 



- 1 



(54) 



On the left hand side, we apply the trivial identity ah — 1 = a(b — 1) + a — 1 



ru+e 



with a - -^y 
for £ —7- 0. Thus 



and b = E 



,e la( T -cc) 



we divide by £ and we take the limit 



0(6-1) r (§ + £)£ 



H T -oJ 



- 1 



r(l) 

^ EMTtJ] 



and: 



0-1 i/ r(} + e) \ 1 / rQ + e )-rQ) 

^ -Lr Q) = -L (-^F) (c s + 2 In 2) = - (c E + 2 ln2) . 



On the right hand side of ( 1541) . we have: 



(sinh(cz)) z 



1 

£ 



exp e In 



(sinh(cz))' 
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hence: 

1 f 00 dz / (sinh(cz)) 2£ 
1J cosh(^) 1 2^ 



£ ^o f x dz I (sinh(c,2)) 2 

cosh(f) n l 2 



uumi^-^ 



oo 



In (2) + 2/ — ^— (ln(sinhM)) 
/ cosh(^) 



which finishes the proof. 



Remark 2.11 a) VFe denote now: 

r°° d 

/o cosh(f) 



Thus: 

cfe /sinh(cz)\ c _^o /°° dz ln(^) 



F(c) - ln(c) ee [^ d ] , In fggMggT) ^o / 

1 ^ l ^ Jo cosh f V c J Jo 



cosh(f) \ c J J cosh(f) 



a 



f°° dz 

F (°) = / ufnz\ ln ( Sillh M) • ( 55 ) 



-0.7832. 

(56) 



b) More generally, we denote: 

FM)s f ;»" l(si " lh(ra)) ' (57) 

and, changing the variables: z = t^u, we obtain: 

'^-S)f=5R ta (-(4-))-3'( e 3)- < 58 > 
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3 The Ornstein-Uhlenbeck case 

3.1 An identity in law for Ornstein-Uhlenbeck processes, 
which is connected to Bougerol's identity 

Consider the complex valued Ornstein-Uhlenbeck (OU) process: 

Z t = z + Z t - \ / Z s ds, (59) 

Jo 

where Z t is a complex valued Brownian motion (BM), zq G C and A > 
and Tc = T°_ cc = inf {t > : \9f\ = c } {Of is the continuous winding 
process associated to Z) denoting the first hitting time of the symmetric 
conic boundary of angle c for Z. It is well known that |32| : 

Z t = e' xt (z + I e Xs dZ : 

= e~ A4 (B at ), (60) 

where, in the second equation, with the help of Dambis-Dubins-Schwarz The- 
orem, (Bf,£ > 0) is a complex valued Brownian motion starting from zq and 

t 2Xt _ 1 

2As - e ' 



a t = / e 2As ds 



2A 

We are interested in the study of the continuous winding process 9 
~hca.{f Q ^§ M -)it > 0. By applying Ito's formula to (160|) . we have: 

dZ s = e- Xs (-\)M as ds + e- Xs d (B Q J . 

We divide by Z s and we obtain: 

dZ s dM as 

-A ds + — - — , 



Z R 



dZ s \ (dl 

1m — -- = 1m 



hence: 



which means that: 

Of = 9l t 
Thus, the following holds: 

24 



Proposition 3.1 Using the previously introduced notation, we have: 

Of = C (61) 

and: 

T c W = l-ln(l + 2AT!; c ), (62) 

where T_ c c is the exit time from a cone of angle c for the complex valued BM 

B. 

Proof of Proposition 13. II We define 

T C W ee T^^mi {t > O:\9f \=c} 

= inf{t>0: \el\=c}. (63) 

Thus, we deduce that «„(» = Tf cc = T_ . However, T B _ CC (the exit time 
from a cone for the BM) has already been studied in the previous chapter 
and we know the explicit formula of its density function (Proposition 12. 81) . 
Thus: 

*?> = a" 1 (< c ) = a- 1 (T*J , (64) 

where a _1 (t) = ^- In (1 + 2At). Consequently: 

T C M = A. i n (i + 2AT! C)C ) , 
and: 

E[T c W]=i-E[ln(l + 2AT e CiC )], (65) 

which finishes the proof. 

□ 

From now on, for simplicity, we shall take z = 1 (but this is really no 
restriction, as the dependency in z Q , which is exhibited in ([7]), is very simple). 
The following Proposition may be considered as an extension of the identity 
in law (ii) in Proposition 12.41 which results from Bougerol's identity. 
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Proposition 3.2 Consider (Z x ,t > 0) and (U x ,t > 0) two independent 
Ornstein- Uhlenbeck processes, the first one complex valued and the second one 
real valued, both starting from a point different from 0, and call T b (U x ) = 
inf {t > : e xt U x = &}. Then, an Ornstein- Uhlenbeck extension of identity 



in law (ii) in Proposition\2.4\ is the following: 



z x (law) 

T (A) (C/ A ) ~~ U a(b), 

where a(x) = argsinh(x). 

Proof of Proposition 13.21 Let us consider a second Ornstein-Uhlenbeck 
process {U x ,t > 0) independent of the first one. Then, taking equation (ToOl) 
for U x , we have: 

e xt U t x = 6 f jx t _ lv (67) 

where (5 t ,t > 0) is a complex valued Brownian motion starting from Zq = 1. 
Thus: 

^ (A) (^ A ) = ^ln(l + 2ATfl. (68) 

Equation (161 p for t = ^ln (l + 2AT fe 5 ), equivalently: a{t) = T^ becomes (we 
suppose that Zq = 1): 

nZ x _ nZ x _ /)B ('??) ^ 

a 



3.2 On the distribution of T_ n r for an Ornstein-Uhlenbeck 
process 

Now we turn to the study of the density function of: 

TW^T!:^inf{t>0:|^|=c}, 
and its first moment. 
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Proposition 3.3 Asymptotically for A large, for z = 1, we have: 

2XE[T^]-\n(2X) X ^E[\n(T e _J], 

and: 

r°° d? 
E [In (TtJ] = 2 ^ cQsh In (sinh (cz)) + In (2) + c E , 

where Ce is Euler's constant. 

For c < |, we have the asymptotic equivalence: 



1 
A 



E [T C (A) ] - E 



sinh ( B. 



A^O 



-E 



sinh [Br, c J) 



Equivalently: 

d_ 
d\ A=0 



A^O 



^ fc (A) ] = ! im . t(^^ (a) ]-^K 0) ]) 



A 



-E 



sinh f B< 



Moreover: 



E 



sinh ( B 



5^7 



dz 



cosh (^) 



(sinh(c2;)) . 



More precisely, for c < ^ : 



E 



[ sinh f 



B 



and asymptotically: 



T 7 „ 



£ 



cos(4c) cos(2c) 



sinh f B. 



~ 5c 



Proof of Proposition 13.31 
A large 



Let us return to equation (1651) . For A —¥ oo, we have: 
B[l?)] = ^B[ln(l + 2AT! C)C )] 



-E 



2A 

In (2A) 1 
2A + 2~V 



i»(2A(r c>c + - 



£ 



In T 



2A 



(69) 
(70) 



(71) 



(72) 



(73) 



(74) 



(75) 
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Thus: 



2XE[T^]-\n(2\)^E[\n(T e _J], 

which is precisely (169|) . Moreover, by the integral representation (}5Tj) for 

E [In (T! cc )] , we deduce (170]). 

A small 

We shall now study the case A — > 0. We have that: 

T C M = i In (1 + 2AT! CiC ) . 

For c < |, from Spitzer ([3]), (at least) the first two positive moments of 
T- cc are finite: E [(T^ cc ) ] < oo, (p = 1,2). We make the elementary 
computation: 



1 /ln(l + 2Ax) 



A 



2A 



x 



1 I 1 



A V.2A J 1 



l+2Az 



' dy_ 

y 



x 



1 r^Xx / i 
y=l+a 1 / /I 



2A 2 



., , , a=2\b n I X b db A^O 2 

1 I da = -2 / rrr > -X 1 . 



f \l + a y i 1 + 2A6 

Consequently, by replacing x = T^ cc , we have: 

r-T 9 



A 



(e[tW]-e[t! c J)=e 



bdb 
1 + 2A6 



We may now use the dominated convergence theorem [7J, since the (db) 
integral is majorized by (X^ cc ) 2 , which is integrable. Thus: 

\{E[T^]-E[T s _J) X -^-E[(T^ c f]. 



Following the proof of Proposition 12.101 Bougerol's identity and the scaling 
property yield: 

E [(sinh (B t )) 2 ] = E [(5 At ) 2 ] = E [A t ( 0l ) 2 ] = E [A t ] E [(5^] 
= E[A t }. 

Thus, for t = TZ„ „, we have: 



E 



An 



E 



sinh f B* 
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Similarly: 



E [(sinh [B t )f] = E [(6 A f] = E [{A t f [5^] = E [{A t f] E [{6^} 

= 3^[(A) 2 ]. 



Thus, for t = TZ C c , we have: 



E 



An 



-E 



( sinh ( B 7 



So, because A T ~t = T_ cc , we deduce (17T1) . In order to prove (1721) . it suffices 
to remark that: 



E pTf] = E [TlJ = E 



A 



T 7 „ 



E 



sinh ( B t j 



On the one hand, by using the density of B T i 



E 



sinh f B T i 



dy 1 

2c cosh(f ; 



(sinh (y)) 1 



i i 



oo 

dz- 
o cosh( 



2 ' 



(sinh(cz)) , 



which is finite if and only if c < |. In order to prove this, it suffices to use 
the standard expressions: sinh(x) = eX ~^ and cosh(a;) = eX+ ^ - . On the 

other hand (note T = TZ CC ), we remark that — B? = B? and [32J, ex.3. 10, 



E [e kBT ] = E 



e 2 



T 



1 



cos(fcc) 



for < k < tt(2c)-\ thus: 



£[(sinh(I? T )) 4 ] 



2' 1 
1 

¥ 
l 



-E 



( e B T _ e -B T y 



E \ e iBT - 4 e 3B T-z T + Q e 2B T -2B T _ 4e B T -3B T + e - 4B T~| 



4 (2E [e iBT ] - 8E [e 2BT ] + 6) 



2 1 
1 



1 



4- 



1 



2 3 Vcos(4c) cos(2c) 
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3 , 



which is precisely (I74p and this is finite if and only if c < ~. Moreover, 
asymptotically for c — ¥ 0, by using the scaling property, we have: 



E 



sinh ( B 



E 



c 4 3E 



sinh ( cB> 



e->0 



c*E 



B 



T' 



(^l,l) : 



5c 4 , 



{Tl t 1) = 5/3 (see [31j; by using the notation of this paper, Table 
for X t = Ci and t = 1). This asymptotics may also 



since E 

3: E [x\] = ^M 

be obtained by (I7^|) by developing cos(4c) and cos(2c) into series up to the 

second order term and keeping the terms of the order c 4 . 



□ 



Remark 3.4 If we slightly modify the above study for the Ornstein- Uhlenbeck 
process by inserting a diffusion coefficient D: 

Z t = z + V2DZ t - A / Z s ds , 

Jo 

we obtain: 



-XI 



-xt 



z + V2D I e Xs dZ t 
o 



(76) 



where in the second equation we used Dambis-Dubins-Schwarz Theorem with 



,2 At 



a t = 2D e 2As ds = D - 

'o -* 



1 



a 



In! l + 4<). 



. -1 __ J_ 
1 ~ 2\'" V ' D 



Thus: 



2\E[T^]-ln(±)^E[ln(T°_J}, 



(77) 
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because: 



E[T(»] = ±-E 



2A 

2A 



2A 2A 



In 1 + — T b 
1 D " 



'■^K + T 



In I Ti v + £- 



Moreover: 

E[ln{T e _ CtC )} = 2ln(z )+E 

= 2\n(z ) + 



In T 



-,0(1) 



dz 



cosh (f ) 



In (sinh (cz)) + In (2) + c E , 



(78) 



n0(D 



where T_ c c denotes the first hitting time of the symmetric conic boundary of 
angle c for a Brownian motion Z starting from 1. 

2 

For X small, we replace 2T _ C c by j^T e c c in the proof of Proposition Iff. 31 (X 
small case) and we have: 

By repeating the previous calculation, we make the elementary computation: 

2 



! /ln(l + fx 
A 



1 rx 



2A 



D 



-x 



% ) bdb x 



2 -/o 1 + A^6 



g) 



:r 



VKe replace x = T_ cc , and by the dominated convergence theorem [1], for 
c < |, we obtain: 

>2 



I ( B [T«] - ^£ (sinh (^ 



A^O 



i(l)^[( r -«) 2 ] 



~, 



E 



sinh f i? 



'r 1 



where E 



sinh f £> 



c^O by ([75p. 



3 V.2£>, 

zs groen fry ( fyffl) . ( f7^| ) and asymptotically, for 
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